The cosmological model with an analytic exit 

from inflation. 



Chervon, S.V. W Zhuravlev, VM.W 
February 5, 2008 



(1) Department of Theoretical and Mathematical Physics, Ulyanovsk State 

University, 

42, Leo Tolstoy, Ulyanovsk, 432700, Russia 
Email: chervon@sv.uven.ru 

(2) Department of Theoretical and Mathematical Physics, Ulyanovsk State 

University, 

42, Leo Tolstoy, Ulyanovsk, 432700, Russia 
Email: zhuravl@sv.uven.ru 

Abstract 

A cosmological model of homogeneous and isotropic spatially fiat 
Universe with gravitating self-interacting scalar field is considered. 
The exact solution, admitting an analytical exit from inflationary 
stage into a radiation era and a matter dominated epoch, is obtained 
by virtue of "fine turning of the potential" method. We found that an 
inflationary stage is supported by decay of higgs bosons in the frame- 
work of the solution obtained. Freidmann's regim is associated with 
adiabatical expantion of the Universe, filled by the matter with special 
equation of state. 

Thus we presented the exact solution which solve the problem of 
transition from an inflationary to a radiation eras or long standing 
'exit' problem. 



PACS-94: 04.20.-q 98.80.-Dr 
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1 Introduction 



In this contribution we investigate general properties of the evolution of a 
homogeneous isotropic universe on the basis of an analysis of the system of 
Einstein equations for a self-interacting scalar field. One of the important re- 
quirements imposed on early universe theories by present-day observational 
data is the existence of an inflationary phase in the evolution of the uni- 
verse with the scale factor increasing at an exponential rate.jl] Thus the 
morden general consept of universe's evolution include the evolution from 
the Planck's epoch through an inflationary stage which should be ended by 
exit to radiation era and then to matter dominated epoch to the present-day 
universe. In a literature we can find a lot of theories which can describe 
all of mentioned epochs, but separately, that is the the physical content of 
the model is changed from one regim to enother. Therefore it is important 
to construct the model of one physical content describing the global evolu- 
tion of the Universe starting from the preinflationary epoch and ending with 
later epochs in the evolution of the universe. The model of this kind should 
have a solution of 'exit problem' from inflation stage to radiation and matter 
dominated stage, mentioned in the original work of A.Guth as a central 
problem of inflationary scenario. In the present work we present the exact 
solution, admitting an analytical exit from inflationary stage into a radia- 
tion era and a matter dominated epoch. The result is obtained by virtue of 
"fine turning of the potential" method and by the analysis of some corre- 
spondence between exact solutions and approximate ones in slow roll regim. 
We found that an inflationary stage is supported by decay of higgs bosons 
in the framework of the solution obtained. Freidmann's regim is associated 
with adiabatical expantion of the Universe, filled by the matter with special 
equation of state. 



2 Survey of exact solutions in inflation 

The inflationary phase of the (isotropic and homogeneous) universe is cus- 
tomarily understood to mean the period £j < t < tf of its early evolution, 
beginning at the time tj and ending at the time tf, when the scale factor R(t) 
increased at an accelerated rate, i.e., 

R = R(H 2 + H)>0 
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(see |I3]|). In this setting the term subinflation is used when H < 0, standard 
or exponential inflation when H = 0, and superinflation in the case H > 0. 

We take a special look at the case in which the law governing the variation 
of the scale factor R(t) has one of the forms 

R(t) = R t m , m>l,R >0, (1) 
R(t) = R e H °\ H > 0, (2) 
R(t) = Ro(U - t) n , n < 0, U > U, (3) 

which establish the general inflation condition R > 0. The scale factor ([I]) 
corresponds to subinflation and is called power-law inflation. The case of 
superinflation (|3|) is unattainable within the scope of self-interacting scalar 
field theory. 

The existence of an inflationary period of expansion of the universe is 
necessary for solving the horizon (homogeneity), flatness, and relic monopole 
problems, which have been discussed in detail, for example, in the original 
work of Guth (see also Linde's book [|TJ). The specific mechanism driving 
inflation in the theory with a self-interacting scalar field (inflaton field) is 
usually identified with the form of the functional dependence of the effective 
self- interaction potential of the scalar field on the field itself: V = V{<p). An 
important factor here in substantiating the existence of such a mechanism is 
the slow-roll regime, i.e., a regime in which the variation of the field [and 
the potential V(<p)] was relatively slow, so that the "kinetic" energy of the 
scalar field can be disregarded: (f> 2 V(<p) ([pE|). 

In the standard model (§) the self-interaction potential is specified as a 
quadratic function of (p: 

v{<p) = n 2 4> 2 + b. 

For this case the field varies linearly with time: oc t, whereas the scale 
factor increases exponentially: R(t) oc e Ht (|IJ and fL5|). 

A new approach to the construction of exact solutions in inflationary 
cosmology (without any reliance on the slow-roll regime) has been developed 
in [12|, [16|] , and |J . The method can be broadly summarized as follows. We 
specify the regime of evolution of the scale factor R = R(t) and determine 
the evolution of the self- interaction potential V = V(<p) in such a way as to 
ensure the specified R = R(t) regime. We then determine the rate of change 
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of the scalar field <fi(t). Next, integrating over t, we find the law of evolution 
of the scalar field = <p(t), which will also be matched with the regime chosen 
for the scale factor. In the final analysis we obtain a parametric dependence 
V = V((p), which then represents the self-interaction potential "fine-tuned" 
to the given exact solution. 



In previous work, |L2] in particular, the self-interaction potentials have 
been plotted as a function of the scale factor R(t) corresponding to the infla- 
tionary regimes (|l|) and @. Also in |12| V(4>) has been obtained in general 
parametric form for certain types of large- 1 asymptotic behavior of R(t) such 
as to ensure regime (|TJ) or (g). 

A somewhat different approach to the construction of exact solutions in 
cosmological inflation models has been presented by Barrow. |T7|] The crux of 
his approach is that the evolution of the scalar field <fi = <fi(t) is specified; the 
next step is to determine the evolution of the scale factor R = R(t) and the 
potential, which depends explicitly on the scalar field <fi. 

We also mention the work of Maartens et al.,Hl8| who have used a method 
similar to fine tuning of the potential, but with a special parameter intro- 
duced in place of time. This approach might well be practical for the solution 
of some problems, but it complicates the process of obtaining and analyzing 
solutions. 

Here we follow the general scheme of the potential fine-tuning method 
and submit two new techniques for constructing classes of exact inflationary 
solutions for a homogeneous and isotropic universe. One technique entails a 
variational formulation of the conditions underlying the slow-roll regime, and 
the second rests on the possibility of specifying the self-interaction potential 
as a function of time. We investigate the general implications of such an 
approach and present new exact solutions. We proposed as well the cosmo- 
logical model based on exact solutions in which the quality solution of the 
global evolution of the universe is found. 

The main attention is consentrated on the solution of the exit problem 
from inflation, because the trasition from radiation- to matter-dominated 
epochs have beed well investigated |^|, [2"T| |. 

Let us remember that the exit from inflationary stage should be accom- 
panied by the set of Freidmann's epoch when the scalar factor has the power 
law evolution R(t) oc t s , where s < 1. Besides s = 1/2 - when radiation- 
dominated and s = 2/3 when matter-dominated epochs are realized. Let us 
mention also that each value of s correponds to certain equation of state in 
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standard theory of Big Bang ||. 

For the models of the early univers's evolution, based on the theory of self- 
interacting scalar field, the inflationary regim is not specific one because it 
can be provided by a wide class of the potentials of self- interaction || [7| . Nev- 
ertheless the condition of transition from inflationary regim to freidmann's 
stages is a specific one and can be realized under very special restrictions. 

We whould like to mention here again the work [18] where the inflationary 
model with an analytical exit to radiation stage and dust matter is proposed. 
As a difference from our approuch, the dependence of the Hubbll's 'constant' 
H = R/R of In R is used in ||18|| ; the e-fold number is chosen as one more 
dynamical variable as well. 



3 Exact solution and self-interaction poten- 
tial for the slow-roll regime 

The Einstein equations for a homogeneous isotropic universe, given an arbi- 
trary form of the self-interaction potential V(<f)) of the scalar field <f)(t), can 
be represented by two equations in the class of Friedmann metrics |12] : 



it, n 1 / . R „R 2 2e \ , 




d 2 lnR 2e , 

+ -f^dt + O , (5) 



dt 2 R 2 

where k is the gravitational constant, A is the cosmo logical constant, and <f>o 
is a constant of integration. This system of equations provides the basis of 
the potential fine-tuning method. As mentioned in the Introduction, once 
the evolution of the scale factor has been specified, this method can be used 
to find the appropriate form of the self-interaction potential to achieve the 
given regime of evolution of the universe. 

The method has been used previously JT2], |J to analyze inflationary 



regimes and to demonstrate the existence of a large number of self-inter- 
action potentials admitting such an evolution. The enormous diversity of 
potentials leading to inflation compels us to look for an additional princi- 
ple that might be used to isolate the potential that actually occurred in the 
early stages of evolution of the universe. Considering the special importance 
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of the existence of the slow-roll regime for an inflationary scenario, we look 
into a nonstandard formulation of the slow-roll regime. We define the slow- 
roll regime with the aid of the variational principle of minimum variation 
of the scalar field with variation of the scale factor R(t) and, as a conse- 
quence, minimum variation of the values of the potential V (</)). On the basis 
of Eq. (|) this condition can be specified by the variational equation 




, , d 2 hiR 2e , „ 
«=/'l/~ 5 p- + ^* = '>. <«) 

where U < t\ < t?, < t/. This condition literally means that the evolution of 
the scale factor must be such that the difference in the values of the field </> in 
any finite time interval [t\, £2] will be the smallest among all other possible 
evolutions. The variation of the values of the potential V(4>) is also the 
smallest in this case. 

If we introduce the notation 

I d 2 hjt 2eV' /2 

F <«>= --3r + jp • 



the Euler-Lagrange equations corresponding to the variational problem [13 
assume the form of a pair of equations for F(t) and R(t): 

d 2 \nR 2e 1 



dt 2 R 2 F 2 
d 2 F 2e „ 

Exact solutions of this system are the simplest to find for the case e = 0, 
which corresponds to a Friedmann flat-space universe. In this case we deduce 
the following from Eqs. 0) and (||): 

d 2 lnR 1 d 2 F 



dt 2 F 2 ' dt 2 
The solution of these equations is given by the functions 

R(t) = k (a t + b ) l/a °e cot , F(t) = a t + b , (9) 
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where ao, 6q, 
(D into Eqs. 
the potential: 



Co, and ko are arbitrary constants. Substituting the solution 
(fD and WM, we obtain the following results for the field and 




- — ln(a t + 60) + 
k ao 



A + 3c; 



6c 



e -2a(0-<Ao) _|_ ^£ e -a(0-<Ao) 

a 



(10) 
(11) 



where a = a^y k/2 . It is interesting to note that the solution fl9D-flllD gen- 
eralizes a solution obtained previously ||19|| using the method of generation of 
new solutions by means of invariant transformations of 0, R, and V, which 
do not alter the equations of the standard inflation model. As should be ex- 
pected, the given solution describes an inflationary regime of the exponential 
or power-law type or a combination thereof. 

The solution (0), ([□]) can be associated with the equation of state of 
matter on the basis of the analogy between self-interacting scalar field theory 
and an ideal fluid. This analogy, in particular, is expressed in the fact that 
by comparing the energy-momentum tensors of these two models one can 
formally calculate the pressure p and the energy density p of an ideal fluid 
in terms of the parameters of the scalar field model according to the rule 

P=U 2 + V{<f>), (12) 



T 



-Tl 



rp2 

~ 1 2 



T">3 

~ J 3 



P = ^ 2 -v{<t>). 



(13) 



Substituting the expressions for the self-interaction potential V(<p) from Eq. ( p]) 
and the field <fi from Eq. (|TUD, we obtain the relations 



P 



1 



A 



3c 2 + JL e -2a(^0) + 

a n ao 



'"0 



p= 

K 



A + 3cq + (— - 2)e 
a 



-2a(<j>-<j>o) 



a(0-</>o) 



6co ( 
a 



-a(<j>—cj>o) 



Eliminating the field <fi from these equations, we obtain an effective equation 
of state of matter in the form 

2 



P = ~P + 



2o| 

K 



" C + \ o( K P ~ A ) 



(14) 
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When the class of solutions (0)— (O) is extrapolated to small and large 
times, i.e., beyond the limits of the inflationary phase U <t\ <t2 <tf, the 
following characteristics are readily established. The given solution always 
begins from a singular state, i.e., for any model parameters (A, a , c , b ) 
there is a time to — —bo/ a o i n the history of the evolution of the scale factor 
when R(to) = 0, after which the scale factor increases at an accelerated rate 
for some time. A natural transition to Friedmann expansion does not take 
place at larger times, and this is a drawback of the analyzed solutions. The 
transition to the Friedmann regime requires that l/al = 2/3 and c = 0. It 
is therefore obvious that without the cosmological constant A in Eq. (|J) its 
role is assumed by Co, and the condition cq = corresponds to the standard 
departure from the inflationary regime. 



Again we emphasize that the resulting analytical solution (p|)-(pTD corre- 
sponds to an exactly tuned potential in the slow-roll regime in the variational 
formulation. 

It is a well-known fact that the pioneering work on inflation (see the 
surveys in [1], [14] , and ]T3] ) treated a truncated system of Einstein and scalar 



field equations, with the slow-roll regime defined in such a way as to eliminate 
the second time derivative of the scalar field <fi from the equations, a device 
justified by the condition 2 <C V((f>). There are no significant constraints on 
the self-interaction potential in this approach. In our case we have obtained 
an exact profile of the potential in the explicit form fllTl). We note once 
again that a solution analogous to (|ll]) has been obtained previously, [|19[ but 
without any mention of the slow-roll regime per se. It is easily verified that 
the condition 2 <C V{4>) is satisfied in a certain time period t{ < At < tf 
for the exact solution (|9|)-(|iT|). 



4 Generation of exact solutions for a given 
potential V = V(t) 

The problem of analyzing the interrelationship between the scenario of evo- 
lution of the universe and the form of the self-interaction potential can be 
treated in time scales greater than the inflationary phase. In this section, 
with a view toward completeness, we discuss the overall evolution of the 
universe. 
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In the standard approach the Einstein equations contain the self-inter- 
action potential of the scalar field in the form of the function V = V(4>). 
In our approach V(<f>) is actually replaced by the function V(t), which must 
be interpreted as an effective potential of material fields. We note that in 
the given case of a homogeneous and isotropic universe the representation 
of the self-interaction potential in the form V = V(t) does not conflict with 
the general variational problem of the derivation of the Einstein equations, 
where significant use is made of the functional dependence V(4>), since each 
function V(t) and <p{t) corresponds one-to-one with a particular function 
V(4>). We call the function V(t) the evolution of the potential energy or 
the potential history, underscoring the departure of our approach from the 
standard approach, in which the dependence V = V(<fi) is fixed. 

We consider the problem of finding all possible types of evolution of the 
universe for a fixed potential history V = V(t) describing the time variation 
of the potential energy of material fields. Stated in the indicated form, the 
problem has exact solutions, whose construction reduces to the integration 
of linear equations in the case e = 0. We write Eq. (^) as an equation for the 
function Z(t) = Z 3 : 

- Z + 3(KV(t) - A)Z - 6eZ 1/3 = 0. (15) 

In the case of a flat-space universe (e = 0) this equation acquires the form 
of an ordinary linear differential equation 

-Z + 3(KV(t)-A)Z = 0, (16) 

which has the same form as the Schrodinger equation for the motion of a 
quantum particle in one-dimensional space with a particle potential energy 
U(t) = 3nV(t) and self-energy E = 3A. Now the function Z(t) assumes the 
role of the particle wave function. We wish to examine this case in more 
detail. 

We assume that one of the solutions of this equation for a fixed depen- 
dence V = V(t) has been found. We denote it by Zi(t), whereupon a second 
linearly independent solution of this equation Z%(£) for a fixed function U(t) 
is easily found, because any two linearly independent solutions of Eq. (|16|) 
are related by the equation, 

Z,Z 2 - Z 2 Z, = W , (17) 



9 



where Wq is a constant. From this result we obtain 



Z 2 = Zx \Qo + Wo J ^ J , 



(18) 



where Qq is a constant of integration. All solutions of Eq. (|T6|) corresponding 
to a specified fixed potential U (t) can be found by varying the constants Q 
and W Q . 

For example, the solution (|^) for R(t) in the problem with a minimally 
varying field leads to the solutions 



R 2 (t) = Ro(t) 



Qo + W J R~\t)dt 



nl/3 



k {a t + b ) 1/a °e cot Q + k 6 W Q f dt(a t + 6 )" 6/a »e- 6c °' 



1/3 



(.19) 



In the case of pure power-law inflation Co = in the solution (§) with the 
minimally varying scalar field we obtain a general class of solutions for V(t) 
fixed in (0): 



4 p) (t) = k (a t + 6 ) 1/a ° [Qo + k 6 W J dt{a t + b ] 
= k (a t + 6 ) 1/a » feo + ^oA(a f + 6 )" 6/ 



-6/ag 



1/3 
1/3 



(20) 



where A = ^/(ciq — 6). 

Note that the solutions ( [Bp and ( 2D ) are new exact solutions for 
metric dependence V = V(0), i.e., = V(t), = 0(t), indicated in the 
solution (|10D (HID- 

The proposed method based on (|I7D carries over without too much diffi- 
culty to the case e ^ 0. It is necessary here to fix the function 

Q{t) = 3(«y(t) - A) - 6eZ- 2/3 

as a time function. The potential V(t) then depends on the form of the 
solution. The case e ^ requires separate investigation. 

The linearity of Eq. (|16D for the function Z(t) permits the formulation of 
a problem in eigenf unctions and eigenvalues, the role of which is taken here 
by the cosmo logical constant, provided that the given equation is augmented 
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with homogeneous initial conditions. For example, we can investigate all 
possible evolutions for a fixed function V(t) such that evolution begins at 
a certain time to from the state R(t) = Z(t) = and returns to the same 
singular state at another time t\ > to- Such scenarios of the evolution of 



the universe can be set in correspondence with oscillating solutions [^1 
The universe originated and disappeared in the time period (t , ti). The 
corresponding problem appears as follows: 

-Z + 3(nV(t) - A)Z = 0, (21) 
Z\ t=t0 = 0, Z\ t=tl = 0. (22) 

This problem has the same form as quantum-mechanical problems for a dis- 
crete spectrum. If the initial conditions are set at times t = —oo and 
t\ = +oo, and the potential energy is a smooth function of time, the uni- 
verse can evolve from other than a singular state and again evolve into a 
nonsingular state in the sense that the energy densities and other physical 
characteristics of matter (field) are finite. Other types of homogeneous initial 
conditions are also possible. The fact that only a finite or denumerable set 
of admissible values of the cosmological constant occurs for each problem of 
this kind can shed light on the issue of the true value of the cosmological 
constant and the physical reasons for it. 



5 Analysis of the evolution of the universe 
for various types of potentials 

The determination of the special characteristics of the origin and evolution of 
various inflationary regimes can be pursued on the basis of the representation 
( |16D in the example of a series of models that are simple from the standpoint 
of the construction of solutions but are interesting from the standpoint of the 
behavior of the potentials. In the present study we have chosen the following 
model potentials in this category: 

3KV(t) = 2t 2 , (23) 

TYl 

3K,V(t) = — , m = const, (24) 

3kVW - -z^ry (25) 
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We investigate the solutions for the potentials (^3|), fl2"4|), and (^) from 
the viewpoint of the possible existence of inflationary regimes and their tran- 
sition to a Friedmann phase of evolution. 

(a) It is readily verified that the potential ( p3[) admits oscillating solutions 
corresponding to the statement of the problem in the form (|2l|) , i.e., laws of 
the evolution of the scale factor such that the universe goes from a singularity 
in the limit t — > —oo to a new singularity in the limit t — > +00. Thus, the 
solutions of Eq. (0), subject to null boundary conditions as t — > ±00, have 
the form 

Z{t) = H t {t)e- t2 l\ 

where Hi(t) are Hermite polynomials. For example, the simplest solution has 
the form Z(t) = exp{— 1 2 /2}. In this case 

R(t) = d exp{-t 2 /6}, =±J^t + <f> , V{<t>) = m) ~ 0o) 2 - 

V 3k 

The value of the cosmological constant ensuring this evolution regime is A = 
1/3 (in appropriate units of measurement for the given problem). For this 
particular solution the condition 

R = (d /3) exp{-t 2 /6}{t 2 /3 - 1} > 

implies the start of inflation at t > and no exit from the inflationary 
regime. 

Other Hermite polynomials correspond to higher absolute values of A = 
1, 5/3, ... . The evolution regimes in this case are such that the universe 
passes through a singular state several times. It is readily verified that the 
inflationary phase for these oscillating solutions begins at t — > —00, where 
the scale factor has a zero minimum, and ends at a time to corresponding to 
an inflection point of the function R = R(t). 

(b) Potentials of the type (|24l) for m >= 1/4 are interesting because in 
the case A = they describe all possible types of power-law evolution of the 
scale factor, including power-law inflation. In fact, the general solution of 
Eq. ( |21~D for arbitrary m > and A = has the form 

Z(t) = dt a + C 2 t p , (26) 
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where a and (3 are distinct solutions of the algebraic equation x 2 — x — m = 0, 
i.e., 

a = - + ^- + m>0, 0=-- 
It is evident from this result that in the case C\ > 0, C2 > the solution 




( |26|) is positive at t > and has one minimum at a point to > 0; after 
passage through this minimum, power-law inflation begins at a time t± > 
to, asymptotically approaching the t a ^ 3 regime. The asymptotic power-law 
inflation regime corresponds to the requirement a > 3, so that m > 6. We 
note that prior to the time t\ the field is imaginary, and the evolution 
regime physically achieved in the self-interacting scalar field model begins 
precisely at t\. To circumvent this problem, it must be required that the 
conditions C\ > and C2 < hold. In this case evolution begins from a 
singular state at a time ts and makes an immediate transition to power-law 
inflation. 

Remarkably, the only regime of power-law evolution of the scale factor 
that does not lead to the potential is Friedmann expansion with the 
ultimately rigid state of matter p — p, for which R(t) oc t 1//3 , so that Z oc t 
and V(t) = —0 for A = 0. The Friedmann expansion regime R(t) oc t 1 / 3 
sets in as t — ► 00 subject to the condition that as the potential decreases, it 
tends to zero more rapidly than 1/t 2 . An example of this kind of behavior is 
afforded by one of the possible functional dependences R = R(t), based on 
evolution of the form 

Z(t) =z + 



t 2 



l + t 

At t > — 1 the scale factor passes through a minimum value z at the time 
t = 0, after which inflation begins. Inflation ends when the inflection point 
is reached, and then as t — > +00, the evolution asymptotically settles into 
the Friedmann regime R(t) oc t 1 ^ 3 . The self- interaction potential in this case 
is 



re 3re(l + t) 2 t 2 + z t + V 
whence it follows that 

2 1 
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For transition to the Friedmann regime with dustlike matter, i.e., when 
p = 0, it is sufficient that V(t) — > 2/t 2 as t — > oo. For any other equation 
of state of the type p = jp, 7 = const, it is sufficient that V(t) — > m/t 2 , 
m > — 1/4 as t — > 00. The case m = —1/4 corresponds to a radiation-domi- 
nated state of matter. 

The potential fl2~4|) also leads to oscillating solutions for m < 0, which 
correspond to A < 0. For example, the general solution of Eq. ([Tj]) for the 
potential (|24]) in the case m = — 2, A < has the form 



Z(t) = A (k - ~) e H + B (-k - e~ fct , jfe = V /Z 3A. (27) 



The set of solutions of Eq. (|27|) includes one that satisfies the condition 
Z(±oo) = 0. It is the solution corresponding to A = and has the form 
Z = A/\t\, A > 0. 

(c) To discern certain characteristics of the limitations imposed on the 
growth rate by the condition for transition to Friedmann expansion, we give 
additional consideration to the history of the potential (|25|) . Like the poten- 
tial (p3|), it increases as t — > 00, but only to a finite value equal to zero. In 
this case the solution for Z(t) with A < has the form 

Z(t) = A(X - A tanh(A t))e A * + B(X + A tanh(A t))e" A *. 

Here A 2 = — 3A > 0. For A = Ao this potential corresponds to an oscillating 
evolution of the form 

Z ^ = 1 » R (f) = cosh^/ 3 (A t), 
cosh(A i) 

where C is an arbitrary constant. This is a unique solution for the given 
potential at A = Ao, corresponding to a unique bound state. As in the 
case of the potential (|2"5|), it describes an inflationary regime in the interval 
(—00, t ), where t is the inflection point of the function R(t). 
In the case A = A = the solution is the function 

R(t) = Ctanh 1/3 A t, 

which describes the progress of evolution from a singular state at the time 
t — to its asymptotic transition to the stationary state R = C = const as 
t — > +00. 
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In the case A > Ao, A > 0, B = we have the solution 

Z{t) = A(X - A tanh(A t))e A *, R{t) = A 1/3 e xt/3 (A - A tanh(A t)) 1/3 , 

which describes evolution without singularities and with the onset of the 
inflationary regime beginning at a time t > 0. 

The solutions corresponding to A > are oscillating solutions with mul- 
tiple passages through the value Z = 0. It is evident that asymptotic Fried- 
mann expansion does not occur for this potential. 

Comparing the solutions obtained for the three potential types, we arrive 
at the following conclusions. First, potentials that increase as t — > oo [po- 
tentials ( ]2"5| ) and (25)] do not allow any transition to the Friedmann regime 



if the growth rate of the potential energy exceeds the rate at which the func- 
tion mt~ 2 with m = —1/4 approaches zero. Transition to the Friedmann 
regime requires that the potential decrease with time according to a power 
law of the type mt~ 2 with m > or that it increase by an analogous law 
with —1/4 < m < 0. Second, all three potential types exhibit the existence 
of the inflationary regime, corroborating the conclusion that this regime is 
not selective with regard to the form of the potential. Altogether these con- 
siderations indicate that the most preferred models are those of the type 
(|24|), which require slight modifications to ensure transition to the required 
Friedmann regime. For example, one such modification is the introduction 
of a weak dependence of m in the equation for the potential (|24[) on the 
temperature T of matter 

V(t,T) = ^Sl (28) 

such that m(T) > near the minimum of R(t) and that m(T) — > 2 in 
the limit t — > oo. For example, the required dependence m(T) can be one 
with jumps characterizing phase transitions in early universe matter. The 
potential (p]) falls in this category in the case Co = 0. Consequently, these 
potentials ( f24|) satisfy the slow-roll principle in the variational formulation 
set forth in Sec. |3|. 

It should also be noted that the simplicity of the potential equation (f28|) 
is preserved only when it is written in the form V = V(t). If the form of 
the function V = V((p, T) is written on the basis of the solutions obtained 
for R(t) and <p(t), this function acquires an intricate form and depends sig- 
nificantly on the constants of integration of Eq. (|I~6D, which, in turn, are 
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determined by the initial conditions for R(t) and 4>{t). This consideration 
demonstrates the preferability of the representation V = V(t) over V = V(<p) 
for analysis of the model dynamics. 



6 The exact and approximate solutions of in- 
flationary models 

The standard cosmological model of SSF, including an inflationary stage, in 
the case of spatially-flat homogeneous and isotropic universe is described by 
the following system of equations 

H 2 = f fy 2 + V{*)) , (29) 

i + 3H<j> = --^-V(<j>) (30) 

where is a scalar field, V(<p) is a self- interaction potential of the scalar field, 

Hit) = j t \nR(t), 

R(t) is a scalar factor, a? is a gravitational constant. The usual process 
of construction of approximate solutions for inflationary regim consist of 
neglecting of the values <ft an d <p 2 in comparision with H<f> and H 2 . This is 
the basis of the slow roll approximation. As a result of such procedure the 
equations above reduced to the following simple form 

H 2 = f (31) 

3Hi = ~V{<l>), (32) 

in which thier integration is a more simple task. For example, follow by the 
review 



IHI , one can obtain 

V 



R{(j>) = R exp{-ae J —d<p} (33) 



For the potential in the trancased equations (pTD-(p2|) and for thiers solutions 
(|33|) the index "" will be used. 
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We prove here that the exact equations (p9|)-(p0]) can be presented in the 
form of approximate ones (]3"ll)-(j3"2l) with the help of an effective potential of 
self-interaction introduced in consideration. We have the formal relation 

i = U(<f>) (34) 

where U(<p) is a function of <fi. The relation (|33j) is true for any dependence 
cf) = (pit). Let us introduce now an effective potential W(<p) following by the 
rule 

W{(j>) = V{4>) + V(0). (35) 

Function W((f>) represents the total energy of the feild as the function of its 
value. Using fl55|) the equations (p9"|)-(|3"0"D can be transfermed to 



H 2 = f W(0), (36) 

3H<j>= -jrW{4>). (37) 

As easily seen, with the accuracy of the exchange of V((j)) on W(4>) equations 

(|36D )([3~7D are equivalent to trancased equations (]3Tf) , (^). A scalar factor 
in the exact model is 

r W 

R{(j>) = tfoexpj-ae J —d<P} (38) 



or 



t 

R(t) = R exp{- J ^W(t)dt} (39) 

to 



The equations of SSF in the form (|3q)-(|37|) are the starting point of the model 
under consideration admitting an exit to the friedmann stage of evolution. 
Equations (|36|),(|3TD lead to the relation 

y/3^UW 1/2 = -W (40) 

which along with ( |33D describe the relation between functions W(4>) and U(4>) 
in a such a way that just only one of the functions above is arbitrary. By 
solving this equation respective to W one can obtain 

w{*) = ^ (/ utfw + y/Wo) 2 > 0, (41) 
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where y/Wo is a real constant of integration. This means that the total energy 
and, consequently, the effective potential W is of nonnegative value in the 
model under consideration. 

In future we will use the following terminology for the potentials V{4>) 
W(4>). The potential of SSF V((f>) presenting in exact equations fl29D -fl30D will 
call true or physical potential The effective potential W(4>) will call effective 
potential of the total energy (for the sake of brivity - the potential of total 
energy) with the aim to differ it from the effective potentials widely used in 
quantum field theories. 

7 The model with effective Higgs potential 

In order to describe an exit from inflation to Freidmann's regim under con- 
dition t — > oo we have to formulate an asimptotical conditions of such exit. 
The sufficient condition following from the observational data is the vanish- 
ing of the energy of an inflaton field when t — > oo. Thus we can require the 
asimptotical condition 

i(0(t)) 2 + V(0)^O, t^oo. (42) 

The last can be decomposed for two another conditions 

0(t)) = U{(t>{t)) - 0, t - oo, (43) 
V{4>{t)) -> 0, t -> oo. (44) 

One of the simplest models satisfing by the requirement ( (43f) is the model 
in which 

= U(<f>(t)) = uoe-^i /i>0. (45) 

For ([|3p we have 

ct>{t) = <Poo- -e""*. (46) 

Therefore 

U(<P)=^ 00 -<p{t))- 
Inserting the last relation to (0), one can find 

W{<P)=(^{<P-<p oo ) 2 -Jv^\ , (47) 
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R{4>) = Rq exp < — ae 



A(0 



4>o 



4A(0 - C 



(48) 



Rq exp 



ae 



exp{- 



33 



Here 



A = /iv3ae, a 
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A 



Vqo is the value of the potential in the limit t —>■ oo, and <fto is the value of the 
field at the starting moment t^. It is easily to check that for the potential 
( f48D the condition (|44l) is valid for the case if = 0. 

From the relation (|35| ) one can find the physical potential 

\ 2 



V{<j>) 



A ( 



(49) 



This is well known potential of 4 -theory with Higgs properties of sponta- 
neous breaking of the symmetry In accordance with this theory the effective 
mass of Higgs boson will be equal to 



M h = \I\^V 00 + ^. 
The minimima of the potential V(<ft) are located in the points 



±2\/A\/r x +^ 



±2^y/3aV 00 + — 



and maximum in the point 0oo (see Fig. 1). 
In the case = we have 



M h = /i, 



L± 



±2\l — . 
3a3 



(50) 



The set of formulas fl4"7|)-(f4"8"D, (pE9]) and (46) gives the exact solution for a 
scale factor of the cosidered model. By inserting in expression for R(4>) the 
dependence = (pit), one can obtain 



R{t) = R exp 



33 



Uq 



e MCrf exp 



3iU 



8/i 2 
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g It isjgllow frorghe relatigs © thgthe condgons © Qd @ eg 
be cons^f red as a p^ssical eqvpvalent of oleay of Higgs bosons with the mass 
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Mh = /i and with the characteristic time of decay « \ Thus in the given 
model it is the decay of Higgs bosons that provide an inflationary temp of 
expantion. It shoud be mentioned here that this model contains the exit from 
inflation but to the asimptotically static Universe. Therefore this model can 
be considered as auxiliary model admitting an important exit property To 
obtain a correct asymptotics when t —>■ oo the model above shoud be modified 
with the account of the following fact: the transition to the Freidmann stage 
is possible if there exists after the moment of Higgs bosons decay a matter 
in the universe in the condition close to thermodynamical equilibrium with 
the equation of state p = 7s. In this case an adiabatical expantion of the 
universe will be of Freidmann type. Using the suggestion above we consider 
the next auxiliary model. 



8 Models with exit to Freidmann 's evolution 

In the work |?j (see also the section 3.) have been analysed the evolution's 
regim's corresponding to the condition 

Vm)) = £f 2 (52) 

where m is the constant, which satisfy to the restriction m > —1/4. In this 
case the scalar factor's evolution is determined by 

R(t) = t^ 3 , (53) 

where 

The equations of states can be obtained from the formula 

(3 — m 1 - 2m + a/1 + 4m 



j3 + m 1 + 2m + y/1 + 4m 



P = 7 £ > 7 

Generally speaking, the solution (|53D can be considered with arbitrary values 
of P > including a power law inflation under j3 > 3. But we restrict 
ourselves by Freidmann's regims when < (3 < 3. The case of power law 
inflation has been considered in the works [§,§], 0- 
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The limiting case m = —1/4 corresponds to equation of state with 7 = 3, 
the case m — - equation of ultrastiff matter with 7 = 1, the case m = 2 

- dust matter with 7 = 0, and the case m = 3/4 - pure radiation with 
7 = 1/3. In general all Freidmann's regim with (3 < 1 correspond to condition 

— 1/4 < m < 6. The model under consideration corresponds to adiabatical 
expantion of the universe filled by the matter with an equation of state close 
to equilibrium one. The special interest is to the radiation era p = | and 
matter dominated epoch p « . 

In the framework of our approach if we have chosen the function V(t) as 
(0) then 

U(t) = V - (54) 

Under this functional dependence U on t the conditions ( f43|) and (|44|) are 
true as well. Following by (|34|) ,(|35D, (|36|) , we obtain the relations 

4>{t) = (pi + z/lnt, U(4>) = z/exp{ (0 — 0i)}, <\>\ = const, 

v 

1 \ 2 
o i/ 2 exp{--(0-0i)} j , 

^z/ 2 exp{-^(0- 0i)}^J - iz/ 2 exp{-^(0-0i)}. 
By inserting in expression for V(4>) the functions (f) = <p(t) we obtain 

V(t) = 

The dependence R = R(t) is defined from (|39|) and for the case under con- 
sideration it has the form 





R(t) = R t^^ 2 exp{ ] I^V 1 t} (55) 

All possible power law scalar factor's evolution corresponding to the depen- 
dence of the potential energy on time (0), can be obtained in the case V\ = 0. 
It is the case that one has an axit on Freidmann's regims ([53]), besides 
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We have to mention, that in the model under consideration, based on 
( |52|) and (54), an inflation is absent, but Freidmann's regim will appear 
asimptotically when t — > oo for any other evolution's regim differing from 
( |52|) and (|54]) by existing of fast decreasing terms. It is just necessary that 
the speed of decreasing is over t~ 2 in fl52| ) and t" 1 in (|54]) , for example, by 
an exponential rate. This point is the basis for the new model generalizing 
the results obtained in sections |7] and |[ 



9 Inflationary model with the exit to the Frei- 
dmann's expantion. 

To construct the model of univers' evolution with inflationary stage admit- 
ting the exit to radiation dominated stage and matter dominated epoch let 
us consider the linear superposition of the exialary models, which have been 
considered in the sections |8|. Namely, let us consider the model corre- 
sponding to the relations 

^ = U{t)=u Q e~^ + ^ (56) 

0(t) = C - ^e-A* + vlnt. (57) 
/i 

In this model the conditions (f43|) are satisfied. Given model has three pa- 
rameters the meaning of which are clear from preveous analysis of the models 
(f4~5|) and (|54]). An inflation regim is possible if the period in the history of 
the evolution exist when 

u e'^ » - (58) 

t 

Evidently, we can choose the parameters fi, v and uq in such a way that the 
condition (|58|) will be satisfied. Let us show this by direct calculations. 

It is difficult to obtain a direct dependence for the function U = U{4>) 
using the relations (|56|) and (|57|). Nevertheless the presentation in the form 
of function U — U(t, 0) is possible. For example, 

U(4>, t) = - + vy\xvt - fi(<f) - (f) c ) 

In the present form we can control the deviation the model (|56D from the 
model ([54]). 
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Using the relations (|40|), (|56|), fl57j ) one can calculate the effective potential 



2 ,. * --frf 



W(t) =l^V + J U\t)dt = e - 2 ^ _ V - + 2n z/ / ^- dt (^9) 



2/x 

to / \ to 

VW = t " f^' - 7 + W ' " + $» 

We will investigate the potentials W((p) and V((f>) as the functions of time, 
because it is impossible to find the explisit dependence on 0. For this purpose 
it is suatable to use the relation (|39|). 

The change of scalar factor R(t) versus the time is shown in the Fig. 3 
and Fig. 4 for the models (|56D, (|57|). Fig. 3 corresponds to the period in 
conventional units before an inflation. Fig. 4 corresponds to the case of 
general evolution of the universe. 

The chois of parameters is special for a demonstration of the main effect 
and equals to v = 6, uo = 30, /i = 0.3, Vq = 0, ae = 1. The change of 
= cf)(t) and V = V(t) versus time are presented in the figures 5 and 6 using 
the parameters above. 

The plot of dependence V = V((p), calculated by plots of dependence 
= (p{t) and V = V(t), is presented in the Fig. 7. 

Thus during the whole evolution the potential is a positive function and 
conseiquently the energy dominant condition is true. 



10 Conclusion 

We have found new classes of exact solutions for a self-consistent system 
of gravitating scalar fields with self-interaction within the framework of the 
cosmology of a homogeneous isotropic universe. We have analyzed aspects 
of the given model that have bearing on the determination of the physical 
conditions restricting the admissible form of the self-interaction potential and 
on the determination of the general form of evolution of the scale factor from 
the standpoint of the existence of an inflationary phase in the evolution of 
the universe, along with other fundamental characteristics of the potential. 

1. We have shown that the definition of the slow-roll regime actively uti- 
lized in inflation is amenable to a variational formulation. In this formulation 
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we have obtained an exact solution for the standard model of inflation with 
a self-interacting scalar field and have given an equation of state of matter 
corresponding to this solution. 

2. We have proposed a method for generating exact solutions with a 
fixed history of the potential V = V(t). On the basis of the method we 
have formulated an eigenfunction/eigenvalue problem for a Friedmann flat 
universe with the role of the eigenvalues taken by the cosmological constant. 
We have shown that in order to analyze the character of the evolution of 
the universe, it is simpler to analyze the model dynamics with the potential 
represented in the form of its history V = V(t) than when it is written as a 
function of the field V = V(<p). 

3. We have analyzed several characteristic types of potential history and 
the corresponding evolutions of the scale factor. Our analysis shows that 
the formal analog of the Einstein equations in the form of the Schrodinger 
equation ([T^) or (|I6D as proposed in this paper can be used to analyze in 
detail the behavior of various physical factors in self-interacting scalar field 
models and to discern a potential selection criterion proceeding from physical 
notions regarding the character of the evolution of the universe in large time 
scales, including inflation (subinflation) as one of the stages, along with the 
conditions for transition of the universe into a Friedmann regime. It follows 
from our analysis that the most realistic model of the history of the potential 
energy is a model of the form ( p8| ) . 

4. We have constructed the model of the evolution of homogeneous and 
isotropic universe posessing an inflationary stage with an analitical exit to 
the radiation dominated epoch and matter dominated era. The model is 
constructed as a superposition of two auxiliary models, the first of them is 
responsible for an inflationary stage, while the second - for the Freidmann's 
regim including matter dominated epoch. 

This work has been carried out with partial financial support from the 
Russian Foundation for Basic Research (Grant No. 98-02-18040). 
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